Ab initio calculation of resonant X-ray scattering in Manganites 
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We study the origin of the resonant x-ray signal in manganites and generalize the resonant cross- 
section to the band structure framework. With ab initio LSDA and LSDA+U calculations we 
determine the resonant x-ray spectrum of LaMnOa. The calculated spectrum and azimuthal angle 
dependence at the Mn K-edge reproduce the measured data without adjustable parameters. The 
intensity of this signal is directly related to the orthorhombicity of the lattice. We also predict a 
resonant x-ray signal at the La L-edge, caused by the tilting of the Mn06 octahedra. This shows 
that the resonant x-ray signal in the hard x-ray regime can be understood in terms of the band 
structure of a material and is sensitive to the fine details of crystal structure. 
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In strongly correlated 3c? transition metal oxides 
(TMO's) orbital degeneracy and orbital order have a 
large effect on magnetic, structural and electronic prop- 
erties. The discovery of a Colossal Magneto-Resistance 
in perovskite manganites [0 has intensified experimen- 
tal and theoretical efforts to understand and quantify 
the repercussions of orbital degeneracy on the physical 
properties of TMO's. An orbitally degenerate ground 
state is unstable and different mechanisms have been pro- 
posed for the lifting of the orbital degeneracy of the 3d 
electrons. One possible mechanism was implicitly pro- 
vided by Goodenough and Kanamori Q, who realized 
that superexchange interactions, which are due to the 
strong Coulomb interactions between electrons, depend 
on the orbitals the spins actually occupy. This implies 
that magnetic and orbital order are directly related to 
each other, and ordering of the orbitals, induced by the 
magnetic exchange, can lift the orbital degeneracy ||. 
On the other hand the lifting of the orbital degeneracy 
can also be induced by lattice distortions: this is the well 
known Jahn- Teller effect j|. In order to infer the ac- 
tual orbital order in TMO's one can therefore make use 
of indirect information coming from the magnetic and 
(Jahn- Teller distorted) crystal structure. For instance in 



), cuprates (KCuF 3 §) 
\2 ] ) the orbital order can be 



vanadates (V 2 3 @, LiV0 2 
and manganites (LaMnOa 
traced in this way. The orbital structure is very impor- 
tant as it strongly affects the low energy and low temper- 
ature behavior of these systems. The description of the 
metal-insulator transition in V2O3 in terms of a spin 1/2 
Hubbard model the magnon dispersion in (doped) 
manganites || and charge-order in half-doped mangan- 
ites are all related to the underlying orbital structure 
of the material. 

X-ray techniques, which could in principle serve as an 
independent direct probe of the orbital order of the 0! 
electrons, are under normal conditions not very sensi- 
tive to the distribution of valence electrons. Neverthe- 



less, it was recently demonstrated on perovskite man- 
ganites that the use of hard x-rays corresponding to the 
Mn fsT-edge, greatly enhances the sensitivity to the va- 
lence electron distribution and is related to the orbital 
order . This type of experiment is now widely used 
to probe orbital and charge order in various TMO's Jl2| . 
The question is, however, where the sensitivity to the 
valence electrons comes from. Essentially three models 
have been proposed. Fabrizio et al. Q suggest that for 
the manganites an Is electron, excited at the K-edge res- 
onance, makes an optically forbidden quadrupole transi- 
tion into the 3d valence band and directly probes the 
empty and occupied Mn e g orbitals. Ishihara et al. Jli[ ] 
argue that the transition is of dipolar character, into the 
Mn Ap states. The sensitivity to the 3d states is then 
caused by the Coulomb repulsion between a Ap and 3d 
electron, which depends on the symmetry of the orbital 
that is occupied. Finally, Elfimov et al. suggest that 
the transition is dipolar, but that the signal is mainly due 
to the in-equivalence of the Ap bands due to the Jahn- 
Teller distortion and therefore not directly sensitive to 
the Mn 3c? states. 

The aim of this paper is to clarify which physical ef- 
fects contribute to the resonant x-ray scattering (RXS) 
signal in manganites. In order to do this, it is necessary 
to make firm theoretical predictions for the RXS signal 
within a chosen model on a specific compound. We adopt 
the model of Elfimov et al. and generalize the resonant 
cross-section to the band structure picture. By means of 
ab initio LSDA and LSDA+U calculations for LaMn03 
we determine the RXS spectrum and show that the cal- 
culated spectrum at the Mn-edge agrees in detail with 
experiment and reproduces the measured data without 
adjustable parameters. This implies that the Jahn- Teller 
distortion is the main cause of the RXS signal. Via a 
tight-binding description of the bandstructure, we show 
that the intensity of the RXS signal is directly related 
to the orthorhombicity of the lattice. We also predict a 
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FIG. 1. Schematic representation of the crystal distortions 
of LaMn03 . Manganese atoms are denoted by the small filled 
circles, the oxygen by the large open circles and the lan- 
thanum by the hatched circles. Left: projection on the ac 
plane of the tilted octahedra. Right: projection on the ab 
plane of the Jahn-Teller distorted octahedra with the occu- 
pied 3x 2 — r 2 and 3y 2 — r 2 orbitals. 

resonant signal at the La L-edge, which was recently ob- 
served in experiment fl(| . The presence of a signal on the 
La-edge is a consequence of the tilting of the MnC>6 octa- 
hedra and, to a lesser extent, of the Jahn-Teller distortion 
of the MnC>6 octahedra in the crystal and is independent 
of the electronic configuration of the Mn 3d levels. This 
shows that the RXS signal is sensitive to the fine details 
of the crystal structure of a material. 

LaMnC>3 is a prototypical orbital ordered compound. 
The Mn 3+ ion is in a high spin d A configuration, where 
three electrons are in the t^g and one electron in the e g 
orbitals. The orthorhombic Pbnm crystal structure can 
be considered as a cubic perovskite with two types of dis- 
tortions: first a tilting of the MnOg octahedra, so that 
the Mn-O-Mn angles become less then 180°, and second 
a cooperative Jahn-Teller distortion (shortening of four 
Mn-0 bonds and elongation of the other two, see Fig. [I]). 
The twofold degeneracy of the e g orbital, occupied by 
one electron, is lifted by the Jahn-Teller distortion. Due 
to these distortions two neighboring Mn atoms are in- 
equivalent: on one Mn-site mainly the 3x 2 — r 2 orbital is 
occupied, while on the neighboring site the occupied or- 
bital has mainly 3y 2 — r 2 character. We show below that 
the in-equivalence of the local environment of neighbor- 
ing Mn-atoms (and also of the La-atoms) gives rise to a 
large signal in a RXS experiment. 

In order to calculate the RXS signal from first princi- 
ples, we first present the formulae for the RXS cross sec- 
tion, that we use later on to obtain the resonant spectra 
and azimuthal angle dependencies from the local density 
calculations. We adopt the notation introduced by Ishi- 
hara and Maekawa |l4| and define the structure factor 

pk-k' = p i(k-k')R, v n \ 

where k (k') is the incident (outgoing) momentum of 
the light, Pj a is the dipole operator with polarization 
a, = x,y, z on site i with lattice vector Rj, |0) (E ) is 



the groundstate wavefunction (energy), \ipq b ) an excited 
state with momentum q, energy E^ and band-index b. To 
be specific we consider the (3,0,0) reflection for the un- 
doped manganite LaMnC>3 in the experimental setup [Tl| ] 
at the Mn if -edge and assume that the incoming light is 
fully a polarized. At the if -edge the relevant dipole op- 
erator is 

Pi,a = l>\ a S i + s lP%a, (2) 

where s, creates a Is core hole on the Mn and pj creates 
an electron in the Mn 4p shell with a — x, y, z symmetry. 
A detailed derivation of the cross-section in the a — ► 7r 
channel, where the polarization state of the photon is 
rotated by 7r/2 on scattering, shows that the intensity is 
given by 

= const |sin0A 2 i(0) - cos0A 3 i(0)| 2 , (3) 

where is the scattering angle, which is fixed in exper- 
iment at approximately 7r/6, and <f> the azimuthal angle 
which describes the rotation around the scattering vec- 
tor, 

A ai(<£) = cos 2 <j>{F xz + F yz )- sin 2 (j)(F zx + F zy ) 

- sin <j> cos 4>{F XX + Fyy + F xy + F yx - 2F ZZ )/V2, (4) 

and 

A 3 i((/>) = cos (f>(F xz - F yz ) 

- sin cp(F xx - Fyy + F xy - F yx )/V2. (5) 

In LaMnC>3 the orbital ordered state consists of two Mn 
sublattices, where on sublattice A (B) presumably mostly 
the 3x 2 — r 2 (3y 2 — r 2 ) orbital is occupied, see Fig@. The 
RXS cross section is exactly sensitive to the difference 
between the two sublattices, as 

— F° a0 - F? a0 . (6) 

i£A ieB 

From this expression it is easy to deduce the physical 
interpretation of RXS cross-section in the band-picture. 
F xx , for instance, is simply proportional to the difference 
of the Mn 4p x projected density of states (DOS) on sub- 
lattice A and the Mn Ap x projected DOS on sublattice B. 
The RXS intensity is proportional to the square of this 
difference. Note that off-diagonal elements of the struc- 
ture factor tensor are non-zero due to the tilting of the 
octahedra. 

The electronic structure is calculated with the LMTO 
ASA LSDA method @ and with the LSDA+U 
method |l9[ , that incorporates a mean- field treatment 
of the strong correlations of the e g electrons. All 
the ab-initio calculations are based on the Stuttgart 
TBLMT047 code. We perform the calculations on a 
set of different (virtual) crystal structures for LaMn03. 



2 



c 



1 




20 30 

E(eV) 

FIG. 2. Resonant x-ray scattering intensity at the Mn-edge 
from LSDA(+U) calculations for different crystal structures. 
JT+T: experimental Jahn- Teller distorted and tilted Pbnm 
crystal structure of LaMn03, JT (T): virtual crystal struc- 
ture with only the Jahn- Teller distortion (tilting of the MnOe 
octahedra). The Fermi- level is chosen as the zero of energy. 
Inset: azimuthal angle dependence of the main edge from the 
JT+T+U calculation normalized to the maximum intensity. 



First of all the RXS signal is calculated for the real lattice 
structure, including the tilting (T) and Jahn- Teller (JT) 
distortion, both with and without the on-site Coulomb 
interaction U. In order to study the effect of tilting and 
Jahn- Teller distortion separately, we do the same calcu- 
lations for LaMn03 in the not Jahn- Teller distorted (but 
tilted) Pr 1 / 2 Sr 1 / 2 Mn03 structure and for a virtual crys- 
tal structure that is Jahn- Teller distorted, but not tilted. 

We first present the results for the Mn K-edge in Fig. ||, 
calculated with Eq.(l-6), where the core hole potential is 
neglected, and the core-hole life-time leads to a broaden- 
ing of the signal with a half width at half maximum of 1.2 
eV (15). In LSDA+U calculations we use U=8 eV |o). 
The spectrum for the full calculation (JT+T+U) shows 



that at 10 eV above the absorption edge there is a strong 
satellite (experimentally 12 eV |ll|,|2l|) and weak satel- 
lites 23 and 27 eV above the edge, that are also found 
experimentally at the same energies . The calculation 
for the JT+T system (U=0) agrees with the result of 
Takahashi et al. p2| and is very similar to the JT+T+U 
result, with a slight shift of spectral weight from the main 
peak to the satellite. Electron correlations, that are re- 
sponsible for the orbital dependent superexchange and 
a possible driving force for orbital order j|, therefore 
make a negligible contribution to the RXS intensity at 
the Mn K-edge, as is suggested by Elfimov et al. |1| 
and in agreement with multiple scattering cluster calcu- 
lations |l7j . The azimuthal angle dependence at the main 
peak shows that the signal follows closely a sin 2 <f> behav- 
ior, in accordance with experiment [O], with a phase 
factor that is slightly shifted by 0.0597T. This indicates 
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FIG. 3. Total resonant x-ray intensity as a function 
of the Jahn- Teller distortion strength calculated within a 
tight-binding framework, normalized to the total intensity at 
the experimental value for the orthorhombicity (defined in the 
text) S = 0.130 (dashed line). 



that the main contribution is due to the term F xx — F yy 
in equation (^), which comes from the difference in p x 
(p y ) projected density of states on the two sublattices 
and is obviously connected to the JT distortion. This as- 
signment is supported by the calculation for the system 
without Jahn- Teller distortion shown in Fig.|^. In the 
T+U system the signal is broad, featureless, an order of 
magnitude weaker and almost entirely due to the tilting. 
From the calculation for the JT+U system we conclude 
that the tilting of the octahedra is responsible for the 
observed satellites, but has only a moderate effect on the 
total intensity of the RXS signal at the Mn if -edge. 

An interesting question is how the RXS signal actually 
depends on the strength of the Jahn- Teller distortion. 
From the LSDA calculation we obtain the parameters for 
a tight binding (TB) description of the electronic struc- 
ture for a lattice without tilting and JT distortion, in ac- 
cordance with previously published results ^3| . The RXS 
intensity for the undistorted lattice and U = is of course 
vanishing because all neighboring Mn atoms are equiva- 
lent. A JT distortion modifies the TB parameters, whose 
distance dependence is given by Harrison |24|] and allows 
the resonant scattering at (3,0,0). In FigjS the total RXS 
intensity is shown as a function of the orthorhombicity. 
The JT distortion gives rise to a long (I) and short (s) 
manganese- oxygen bondlength in the ab plane and we de- 
fine the orthorhombicity as 6 = (I — s)/(l + s) so that for 
the cubic perovskite 5 = 0. The experimental value for 
the orthorhombicity at low temperature is 8 = 0.130 ||. 
The resonant intensity smoothly depends on S and is ap- 
proximately proportional to 5 3 ^ 2 . This supports the sug- 
gestion that the temperature dependence of the RXS is 
closely related to the temperature dependence of the lat- 
tice parameters pl[ . 

Finally we present the results for the RXS signal at 
the La L-edge. Fig. [l] shows that due to the lattice dis- 
tortions also neighboring La-atoms have locally different 
surroundings. We calculate the RXS intensity for the 
(3,0,0) reflection in the a — > it channel, for different (vir- 
tual) lattice structures, as we did for the Mn-edge. We 
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FIG. 4. Resonant x-ray scattering intensity at the La-edge 
from LSDA(+U) calculations for different crystal structures. 
JT+T: experimental Jahn- Teller distorted and tilted crys- 
tal structure, JT (T): virtual crystal structure with only the 
Jahn- Teller distortion (tilting of the MnOe octahedra). 

find a signal at the La-edge, shown in Fig. ||, that is 
mainly due to the tilting. From the JT+U calculation 
we see that the combined contribution of the other two 
symmetry breaking effects is small. For the crystal struc- 
ture that only has the tilting distortion we find a signal 
that is comparable to the JT+T+U calculation both in 
intensity and shape. This result is not unexpected as the 
tilting of the MnOg octahedra modifies the local environ- 
ment of La atoms: the apical oxygen atoms are displaced 
towards the La atoms on one sublattice and away from 
them on the other sublattice so that neighboring lan- 
thanum atoms have inequivalent oxygen surroundings. 

We determined, in conclusion, with ab initio LSDA 
and LSDA+U calculations the resonant x-ray spectrum 
(RXS) of LaMnC>3. The calculated spectrum and az- 
imuthal angle dependence at the Mn-edge reproduce the 
measured data without adjustable parameters. Spectra 
for different (virtual) crystal structures show that the 
RXS intensity at the Mn if -edge is mainly due to the 
Jahn- Teller distortion jL5| and that the tilting of the 
MnC>6 octahedra causes the strong satellite structure in 
the spectrum. We find that the intensity of the reso- 
nant x-ray signal at the Mn-edge is directly related to 
the orthorhombicity of the lattice. We also predict a 
RXS signal at the La L-edge, which is a consequence of 
the tilting of the MnOe octahedra in the crystal. This 
shows that the resonant signal at different edges in the 
hard x-ray regime can be understood and predicted by 
calculations in the band-structure framework and is sen- 
sitive to the fine details of lattice distortions, but does 
not directly probe the orbital order of the 3d states. Or- 
bital order, of large enough periodicity to satisfy Bragg 
conditions, might however be visible in resonant experi- 
ments with soft x-rays, where direct transitions into the 
3d states are allowed. 

We thank I. Elfimov, G. Sawatzky and P.J. Kelly for 
useful discussions. 
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